High-pressure lubricity at the meso- and nanoscale 
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The increase of sliding friction upon increasing load is a classic in the macroscopic world. Here 
we discuss the possibility that friction rise might sometimes turn into a drop when, at the mesoscale 
and nanoscale, a confined lubricant film separating crystalline sliders undergoes strong layering and 
solidification. Under pressure, transitions from N — > N — 1 layers may imply a change of lateral 
periodicity of the crystallized lubricant sufficient to alter the matching of crystal structures, influ- 
encing the ensuing friction jump. A pressure-induced friction drop may occur as the shear gradient 
maximum switches from the lubricant middle, marked by strong stick-slip with or without shear 
melting, to the crystalline slider-lubricant interface, characterized by smooth superlubric sliding. We 
present high pressure sliding simulations to display examples of frictional drops, suggesting their 
possible relevance to the local behavior in boundary lubrication. 

PACS numbers: 68.35.Af,68.08.De,62.10.+s,62.20.Qp 



INTRODUCTION 

Confined lubricants under shear and high pressure dis- 
play, both experimentally [IHS] and theoretically [IHS], 
intriguing nano- and meso-scale tribological phenomena. 
The intervening lubricant film between two sliding solid 
surfaces generally changes from liquid (with hydrody- 
namic lubrication), to solid or nearly solid at high pres- 
sure when the film is only a few monolayers thick. Both 
experiments and simulations find that in this regime the 
film develops a solid-like layered structure, supporting 
static friction, and a strong stick-slip frictional behavior. 
The stick-slip is often believed to be associated either 
with a lubricant melting-freezing mechanism, where the 
film disorders and melts at slip and solidifies during the 
stick phase (see [8] and references therein) , or else with a 
layer-on-layer shear (a shear band) within the film bulk 
or at the film- wall interface [TT3T412J . Moreover, as lubri- 
cant layers get successively squeezed out under increas- 
ing pressure, each N — > N — 1 ( "relayering" ) transition 
generally corresponds to an upward friction jump [21 H]. 
Virtually everywhere in sliding friction, including the 
smoothest boundary lubricated systems, friction grows 
with load. The question which we raise here is whether, 
and under which circumstances, that friction jump could 
be downward rather than upward at relayering. The rela- 
tive surface structure properties at the interface between 
sliders and confined film may play the key role. Under 
pressure the crystal-like lubricant ordering perpendicular 
to the interface will typically be accompanied by ordering 
parallel to the interface [13 H3] . However, little is known 
about possible two-dimensional (2D) parallel crystalline 
order, and even less about its changes upon high pressure 
relayering. There is some evidence, for example in molec- 
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FIG. 1: (Color online) Sketch of the three-dimensional fric- 
tional model (a) and snapshots of three different sliding states: 
(b) 8-layer crystalline lubricant (sticking), (c) 8-layer melted 
lubricant (slipping), (d) 7-layer (relayered) crystalline lubri- 
cant, with a superlubric film- wall interface. The grey-red col- 
ored scale highlights particle kinetic energies, from low to high 
values, (e) Bulk lubricant equation of state (sketch) with 
jumps indicating pressure-induced relayering transitions un- 
der confinement. 



ular squeezout simulations [15] indicating crystalline par- 
allel ordering of the narrowly confined lubricant, and its 
dependence on the load. That detail is important for 
the case of crystalline sliders, for in general the frictional 
forces will depend on the relative geometry and com- 
mensurability at the slider-lubricant interfaces, in a way 
which is not immediately predictable, in particular not 
automatically monotonic with pressure. To explore the 
variety of possibilities, we decided to abandon the open 
geometry where lubricant squeezout takes place, surely 
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much faster with sliding than without [TH HZ] > an d we 
simulated instead lubricated sliding friction under sealed 
conditions, with constant particle number and no squee- 
zout. Sealed conditions might be locally realized in some 
real cases, but we use them here essentially as a theo- 
retical device. Our results show that in agreement with 
pressure-induced N — > N — 1 relayering of the solidified 
confined lubricant film, the parallel lubricant periodicity 
also changes; and as a result the lubricated sliding friction 
may indeed switch from stick-slip, accompanied by bulk 
shear with or without melting-freezing, to smooth sliding 
regimes typical of incommensurate interfaces exhibiting 
superlubric dry friction. In that case the friction coeffi- 
cient correspondingly shows a strong, if counterintuitive, 
non-monotonic behavior as a function of increasing load. 

MODEL 

Our model system comprises two parallel, rigid, pe- 
riodically corrugated plates confining between them a 
small number of lubricant layers, with planar periodic 
boundary conditions (PBC), with n, n t and ny, defining 
the total numbers of particles in the lubricant, top, and 
bottom plates respectively. We will present molecular dy- 
namics (MD) simulation results and a detailed statistical 
data analysis mostly for 2D systems, where the two plates 
consist of onc-dimcnsional line-like confining substrates 
with PBC applied only along the sliding x-direction. We 
also simulated, not as extensively, more realistic (and 
computationally more expensive) three-dimensional (3D) 
systems such as that depicted in Fig. [T|a) observing 
an essentially identical behavior (details in supplemen- 
tal material). The bottom plate is assumed to be rigid 
and immobile while the top sliding plate, whose cen- 
ter of mass coordinate is defined in two dimensions by 
R top = (Xtop, i?top), is connected to an external spring 
of stiffness k driven horizontally along x at constant ve- 
locity vq. A downward vertical load force L is applied to 
the top plate center of mass, Fn — L/n t being the load 
per top slider particle. The corresponding friction force 
(per top particle) is directly measured by the instanta- 
neous spring elongation F = k(X top — u t)/n t . The wall- 
lubricant (WL) and the lubricant-lubricant (LL) interac- 
tions are modeled by Lennard- Jones (LJ) potentials Ull 
and Uwl, choosing an amplitude ratio cll/^wl smaller 
than 1 (typically we used 0.2), so that the lubricant wets 
the substrate, favoring epitaxial ordering in the solid 
phase, such as would generally be the case for example 
for an octamethylcyclotetrasolixane (OMCTS) lubricant 
film between crystalline sliders. We assume aww = 1 as 
our reference length and define <twl = {&ww + cfll)/^- 
Acting on the LJ interparticle distance ratio oll/^wl 
we can explore the effect of wall-lubricant interface com- 
mensurability changes on the tribological response of the 
confined system. The equations of motion for the top 



plate and the confined lubricant particles are 
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where r^, rj, and are the lubricant, top and bottom 
plate particle coordinates, respectively, m is the lubri- 
cant particle mass, M = n t m that of the top plate, 
and i and k are unit vectors along x and z respectively. 
During simulated sliding, the shear-induced Joule heat 
is removed by means of a viscous damping proportional 
to the relative velocity between the top plate and the 
lubricant, r\ being the friction coefficient. Importantly 
here, as detailed later, notwithstanding this rather arbi- 
trary choice of dissipation, picking another thermostating 
procedure only influences details, but it does not change 
the phenomenology essence. The temperature T is con- 
trolled by an ordinary Langevin thermostat, with a ran- 
dom force obeying the fluctuation-dissipation theorem, 
i.e. (f[ an (t)} = and (f™"(t)f™"(i')) = ^]K B T8 l2 8{t - 
if), Kb being the Boltzmann constant. For the top 
plate F ran = -£r=i f r" so that (F ran (t)) = and 
(F ran (t)F ran (t')) = nAr)K B TS{t - t'). 

LJ units are used throughout the paper. The equa- 
tions of motion are integrated using a modified Vclocity- 
Verlet algorithm with a sufficiently small time step (At — 
0.005). Sliding simulations are performed at temperature 
KbT = 0.1. The external driving velocity is vq = 0.1, 
the spring constant k/n t — 0.1, the damping coefficient 
i] = 0.2 (underdamped regime), and the lubricant par- 
ticle mass m = 1. Typical system sizes comprise from 
hundreds up to a few thousands particles. Statistical val- 
ues of the considered physical quantities are obtained by 
averaging over sufficiently long time intervals in steady 
state regimes. 

RESULTS 
Load-induced drop of friction 

We start off with a film of N = 8 solid lubricant layers, 
zero load, and zero sliding speed. Choosing our temper- 
ature above the melting point, the lubricant film first 
melts. From this point we increase by steps the load Fn, 
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FIG. 2: (Color online) (a): average friction force versus adi- 
abatically increasing (circles) and decreasing (squares) ver- 
tical load, for the lubricant-slider 2D commensurate system 
{crLt/crwL = 1)- A hysteresis loop appears at the relayering 
transitions of the lubricant. The numbers 1,2 and 3 refer to 
three different values of load highlighted in Fig. [3] (b) : load 
dependence of the lubricant film thickness, measured by the 
top wall vertical coordinate, {Z t0 p)- (c)-(d): average friction 
force, same as (a) now for 2D lubricant-slider incommensu- 
rate system ((Tll/cwl ~ 1.11475). The 8 — > 7 relayering 
transition here takes place at smaller load, still leading to a 
friction drop. 
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FIG. 3: (Color online) Upper panels: friction force as a func- 
tion of time for the three different increasing loads (Fjv = 1.6, 
3.0 and 3.5) highlighted in Fig. [5] (a), for the 2D commensu- 
rate system (ull/o'wl = 1)- At Fm = 3.5 (panel 3) the 
lubricant film has already undergone the relayering transition 
from 8 to 7 layers. Lower panels: corresponding top plate 
and layer-averaged lubricant velocities, color-coded along z 
according to the inset sketch. At the relayering transition 
sliding switches from the lubricant middle (panels 1,2), shear 
melting with an almost laminar flow, to the hard solid wall- 
lubricant interface (panel 3). 



thus exerting a pressure P; at first the lubricant recrys- 
tallizes with N = 8 layers. Increasing further the pres- 
sure, the compressed confined lubricant roughly follows 
(with a kinetics that is fast here given the small system 
size) the Lennard- Jones equation of state, as sketched in 
Fig. [ije). Specifically, the averaged perpendicular inter- 
particle spacing diminishes gradually, as highlighted by 
the corresponding gentle decrease of Z top in Fig. |2fb), 
up to the occurrence of a relayering transition. As indi- 
cated, the finite system size determines a series of crit- 
ical pressures Pjv-viv-i, Pzv-i-»iV-2 5 etc, where relayer- 
ing takes place, first from 8 to 7 layers, then from 7 to 
6, etc, corresponding to significant downward jumps of 
interparticle spacing and to a sudden change of the mu- 
tual parallel commensurability between lubricant film and 
sliders. The interparticle space shrinking taking place in 
all directions, at each relayering the mutual commensura- 
bility between lubricant film and sliders undergoes a sud- 
den change. At this static level, our choice ° ww = 1 and 

. . . , <TLL 

the relatively strong plate-lubricant interaction, make the 

plate-lubricant interface initially commensurate in the 



low pressure state. That interface becomes incommensu- 
rate as relayering jumps take place under pressure. When 
successively the top plate driving speed vq is raised above 
zero, and sliding takes place, the frictional force, and in 
fact the overall sliding habit, displays a strong depen- 
dence upon pressure. While friction naturally increases 
with load, there are large jumps at each relayering, where 
friction now drops, running against conventional wisdom. 

Figures [2ja) and (b) show the non-monotonic pres- 
sure behavior of the average friction force (P) , and the 
corresponding average effective film thickness (top wall 
vertical coordinate) (Z top ), for adiabatically increasing 
(circles) and decreasing (squares) load P/v, and an ini- 
tially commensurate sliding interface (<Jll/o~wl = 1)- 
At small load, with N = 8 solid lubricant layers, the 
sliding is smooth and friction is relatively low. With 
increasing load, the smooth sliding regime is gradually 
(although not uneventfully, as will be described later) 
replaced by stick-slip. Here friction is much higher, cor- 
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responding (for this choice of the interaction parameters) 
to shear induced melting-freezing [SHE] originating inside 
the lubricant, as visually highlighted in snapshots (b) 
and (c) of Fig. [I] At a critical load Ps->-7> however, the 
lubricant, unable to support the excessive pressure, un- 
dergoes relayering accompanied by an abrupt, and tribo- 
logically crucial, change of parallel lubricant crystalline 
periodicity accompanying the drop of perpendicular in- 
terlubricant spacing. In the denser relayered state, two 
things happen. First, the shear induced melting is now 
more difficult. Second, and most important, the inter- 
faces change their mutual commensurability. Thus the 
relayering causes the shear gradient maximum to switch 
from the bulk lubricant middle to the plate-lubricant 
solid interface (Fig. [I] snapshot (d) and movies in the 
Supplemental Material) whereby stick-slip is replaced by 
"superlubric" [HI HH] smooth sliding, with a dramatic 
friction decrease. Upon further increase of load, the stick- 
slip dynamics returns, this time without lubricant melt- 
ing, until Pt^q and the next relayering take place. Cy- 
cling adiabatically the pressure up and down opens up 
hysteretic frictional cycles close to relayering transitions 
(Fig.^a) and (b)). 

Robustness of friction drop 

This observed frictional drop phenomenon might ap- 
pear at first strictly determined by our special choice of 
parameters, and thus not robust; we found that it is not 
so. 

A variation of the LJ parameter ratio <Jll/&wl makes 
the initial low pressure lubricant-wall interface incom- 
mensurate, yet still leading to clear non-monotonic fric- 
tion and downward jumps, as shown in Figure [2] (c) and 
(d). 

The same qualitative result is also obtained changing 
the total number of confined lubricant particles, and by 
moving to a more realistic 3D model. In Fig. |4j the blue 
data points show an almost laminar lubricant flow of the 
confined film corresponding to a high dissipative stick- 
skip regime of motion. The red points highlight the slid- 
ing switch, after pressure-induced relayering, from the 
bulk lubricant middle to the incommensurate top wall- 
film interface, realizing an almost frictionless (superlu- 
bric) dynamical regime (movies in Supplemental Mate- 
rial) . 

Nonequilibrium molecular-dynamics simulations in 
sliding friction, are also hampered by arbitrariness and 
uncertainties in the way Joule heat is removed; in order 
to attain a frictional steady state, a realistic energy dis- 
sipation is generally impossible to simulate reliably, due 
to size and time limitations (2Q1EI]- However, as shown 
by the average friction force trend in Fig. [5j increasing 
Langevin damping or else adopting a different thermo- 
stat scheme (a damping exponentially decaying with dis- 
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FIG. 4: (Color online) 3D lubricant-wall system geometry: 
(a) average friction force versus adiabatically increasing ver- 
tical load, showing a large frictional drop at the 8 — > 7 relay- 
ering transition; (b) z-resolved profile (layer-by-layer) of the 
lubricant velocity x-component. 
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FIG. 5: (Color online) Average friction force versus adia- 
batically increasing vertical load for two distinct dissipation 
schemes: (circles) a constant and heavy damping (rj = 2.0, 
i.e., 10 times larger than that used for the results reported in 
Fig. [2| on all lubricant particles; (squares) a damping expo- 
nentially decaying with distance away from both substrates 
(see Ref. 0). 



tance away from both sliders [5]) only influences details 
here, without changing the essence: the pressure-induced 
switches of the maximal shear zone from the lubricant 
middle to the slider-lubricant interface, entailing a fric- 
tional drop, is robust and occurs with all thermostating 
procedures tried. 

Another relevant point strengthening the results' gen- 
erality deals with the unimportance of whether dissipa- 
tive stick-slip shearing through the middle of the solid lu- 
bricant film will be accompanied by melting-freezing, or 
will just constitute an interlayer shear band. Experimen- 
tally (e.g., in surface force apparatus measurements), one 



cannot directly observe the film dynamics and both sce- 
narios easily achievable in MD simulations (depending on 
the model parameters, geometry and driving conditions) 
represent interesting perspectives. However the physics 
which we describe is not bound to one or the other. The 
switch of shear from center to interface which we pro- 
pose and demonstrate in a particularly simple model, 
describes a more general phenomenon, that occurs in 
our simulations, where forces are rather generic, in ei- 
ther cases. 

The time dependent dynamical frictional forces in 
Fig. [3] show, at loads Fn marked as 1, 2, 3 in Fig.[2ja), an 
increasingly strong stick-slip followed by a drop back to 
smooth sliding corresponding to a harder solid lubricant 
with a different parallel commensurability at the relayer- 
ing transition. As highlighted by the color-coded lubri- 
cant speed along z (lower panels 1,2, Fig. [3]), the shear 
flow profile is almost laminar inside the film, starting off 
as a shear band at the center of the film. At relayering 
(e.g., 8 — > 7 layers), bulk shear disappears, and the speed 
gradient switches entirely to the (now incommensurate, 
and superlubric) slider-lubricant interface (lower panel 3, 
Fig. [3]). Upon further increase of load, superlubric slid- 
ing persists only up to a threshold pressure. Here, despite 
incommensurability, static friction and stick-slip sliding 
reappear, with friction again rising with pressure. Un- 
like low pressures, high pressure melting-freezing is sup- 
pressed, and stick-slip is ruled by inertia [5]. At each 
successive relayering transition (eg. from 7 — > 6 lay- 
ers), friction may, in principle, jumps down (as shown in 
Fig. [5Ja) and (d)) or up depending, respectively, on the 
realization of a more favorable or unfavorable interface 
incommensurability due to lubricant reordering parallel 
to the confining surfaces. 

For the choice of the interaction parameters in Fig. [2] 
when melting takes place, the lubricant density decreases. 
We do observe that density drop in the form of a small ex- 
pansion of the film thickness in mid-lubricant slips (that 
lead to shear- induced melting), clearly resulting in cor- 
responding vertical jumps of the top confining substrate. 
In contrast, the stick-slip dynamics in the presence of 
a damping exponentially decaying with distance away 
from both substrates (square-point curve in Fig. [5]) does 
not display melting-freezing induced by sliding, with the 
top wall vertical Z top -coordinate just fluctuating weakly 
around a mean value and without any observable density 
decrease in mid-lubricant slips (panel (c) in Figs. [6] and 
0. 



Superlubric sliding under pressure 

This pressure induced switch from superlubricity to 
static friction with stick-slip sliding deserves special at- 
tention, because it takes place continuously and within 
the same incommensurate state of the lubricant. To un- 
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FIG. 6: Stick-slip dynamics in the presence of a damping 
exponentially decaying with distance away from both sub- 
strates, at a vertical load Fn = 2.5 just before the 8 to 7 
relayering transition: (a) frictional force pattern; (b) top wall 
horizontal X-coordinate; (c) top wall vertical Z-coordinate. 
Here we do not observe a density decrease in mid-lubricant 
slips, that now leave the film mostly crystalline and do not 
lead to shear-induced melting. 
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FIG. 7: Stick-slip dynamics as in Fig. [6] but at a vertical 
load Fn = 25 approaching the 7 to 6 relayering transition. 
Here, apart from the substrate-distance dependent dissipation 
scheme, the shear induced lubricant melting during slips is 
hampered by the strong confining pressure. 



derstand that, we calculated the static friction as follows. 
Starting from a set of configurations of the film whose 
loads are above the 8 to 7 relayering transition, we first 
removed the external driving and equilibrated the sys- 
tem under the action of increasing loads. For each con- 
figuration, we evaluated the static friction threshold by 
driving the top wall external spring at a very small veloc- 
ity (i.e., 100 times smaller than that previously adopted) 
and measuring the height of the peak preceding the first 
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FIG. 8: (Color online) Time-evolution of the top plate coor- 
dinate X top for 2D incommensurate sliding with cjll /owl ~ 
1.11475) (same as in Fig. |2jc)), for increasing load just after 
the occurrence of the 8 — > 7 relayering transition. Curves are 
offset vertically for clarity. The inset highlights the Aubry- 
like load-induced onset of the static friction for the top plate, 
marking the transition from superlubric sliding to stick-slip. 



slip. Figure [5] shows a clear onset of static friction near 
Fff = 4.7, taking place without any change in the in- 
commensurate slider-lubricant registry. This represents 
an example of the well-known "Aubry transition" 22J in 
one-dimensional systems, where the increasing interac- 
tion between a harmonic sliding chain and an incommen- 
surate periodic potential gives rise to dynamic pinning 
with the onset of static friction. The sliding trajectories 
(Fig. [8| correspondingly switch from smooth sliding to 
stick-slip of strength increasing with load. 



CONCLUSIONS 

In summary we have explored, by simulating lubri- 
cated sliding under sealed conditions, the possibility that 
frictional jumps associated with N —¥ N — 1 relayering 
transitions might be downward rather than upward as 
commonly observed. We find that downward jumps may 
occur due to two factors: (i) high friction with internal 
shear, easier at low pressure and low densitycould dis- 
appear after relayering, when density is higher; (ii) the 
change of interface commensurability entails possible su- 
perlubricity at higher pressure. In experiments so far, 
the friction drop just described is likely removed by lu- 
bricant squeezout, and routes should be considered to 
overcome that. Despite a clear difficulty in realizing 
sealed sliding configurations in ordinary materials and 
sliding geometries, these two elements stand as results 
of general importance, given the relevance of any ele- 
ments that could control, and especially decrease, fric- 
tion. Under the more common open boundary condi- 
tions of boundary lubrication, where relayering also takes 
place, one still could, depending on conditions and ma- 
terials, realize one or both points above. In realistic in- 
homogeneous conditions, moreover, lubricant "puddles" 
might become sealed under pressure [23], and their sud- 
den yielding could be important for the local sliding dy- 
namics. Finally, future developments of artificial sliding 
colloidal systems [23J [2B] might permit soon the realiza- 
tion of sealed sliding, and the verification of our results, 
including, in addition to the downward jumps, pressure 
induced Aubry transitions. 

This work is part of the Swiss National Science Founda- 
tion SINERGIA Project CRSII2 136287\1. The authors 
gratefully acknowledge D. Vanossi for computational re- 
sources; and N. Manini and O.M. Braun for discussions. 



Possible experimental realization 

First let us consider estimates for OMCTS, a reason- 
ably spherical molecule, as a test case lubricant. Using 
LJ parameters for OMCTS as provided by H. Matsubara 
et al. [24] . the critical pressure for the 8 — ?> 7 relayering 
is about 6 — 7 Kbar. That seems too high compared with 
ordinary squeezing loads in surface force apparatus, al- 
though local trapping of sealed lubricant "puddles" might 
still take place at inhomogeneous interfaces [2"5] . 

Mesoscopic systems appear more promising. For ex- 
ample, colloids can be confined in layers, crystallized due 
to their soft screened Colomb repulsion, and also driven 
to slide against laser generated periodic potentials [2"5T - 
|2"7] . They might therefore represent a good trial case for 
a study relayering of high pressure frictional drops. 



[1] R. G. Horn and J. N. Israelachvili, J. Chem. Phys. 75, 
1400 (1981). 

[2] H. Yoshizawa and J. N. Israelachvili, J. Phys. Chem. 97, 
11300 (1993). 

[3] J. Klein and E. Kumacheva, Science 269, 816 (1995). 
[4] J. P. Gao, W. D. Luedtke, and U. Landman, Phys. Rev. 

Let. 79, 705 (1997). 
[5] J. P. Gao, W. D. Luedtke, and U. Landman, J. Phys. 

Chem. B 102, 5033 (1998). 
[6] P. A. Thompson and M. O. Robbins, Science 250, 792 

(1990). 

[7] B.N.J. Persson, Sliding Friction, Physical Properties and 

Applications (Springer, Berlin, 2000). 
[8] O. M. Braun and A. G. Naumovets, Surf. Sci. Rep. 60, 

79 (2006). 

[9] A. Vanossi, N. Manini, M. Urbakh, S. Zapperi, 
and E. Tosatti(2011), Rev. Mod. Phys., in press; 
arXiv: 1112. 3234v3 [cond-mat . mes-hall] . 
[10] Y. Lei and Y. Leng, Phys. Rev. Lett. 107, 147801 (2011). 



7 



[11] J. Klein, Phys. Rev. Lett. 98, 056101 (2007). 

[12] O. M. Braun, Tribol. Lett. 39, 283 (2010). 

[13] See, e.g., P. A. Thompson, M. O. Robbins, and G. S. 

Grest, Isr. J. Chem. 35, 93 (1995), and references therein. 
[14] K. G. Ayappa and R. K. Mishra, J. Phys. Chem. B 111, 

14299 (2007). 

[15] U. Tartaglino, I. M. Sivebaek, B.N.J. Persson, and 
E. Tosatti, J. Chem. Phys. 125, 014704 (2006). 

[16] B.N.J. Persson and E. Tosatti, Phys. Rev. B 50, 5590 
(1994). 

[17] B.N.J. Persson and F. Mugele, J. Phys.: Condensed Mat- 
ter 16, R295 (2004). 

[18] K. Shinjo and M. Hirano, Surf. Sci. 283, 473 (1993). 

[19] M. Dienwiebel, G. S. Verhoeven, N. Pradeep, 
J. W. M. Frenken, J. A. Heimberg, and H. W. Zandber- 
gen, Phys. Rev. Lett. 92, 126101 (2004). 



[20] A. Benassi, A. Vanossi, G.E. Santoro, and E. Tosatti, 

Phys. Rev. B 82, 081401R (2010) 
[21] A. Benassi, A. Vanossi, G.E. Santoro, and E. Tosatti, 

Tribol. Lett. 48, 41 (2012). 
[22] M. Peyrard and S. Aubry, J. Phys. C: Solid State Phys. 

16, 1593 (1983). 
[23] S. J. O'Shea, N. N. Gosvami, L. T. W. Lim, and W. Hof- 

bauer, Jpn. J. Appl. Phys. 49, 08LA01 (2010). 
[24] H. Matsubara, F. Pichierri, and K. Kurihara, J. Chem. 

Theory Comput. 6, 1334 (2010). 
[25] T. Bohlein, J. Mikheal, and C. Bechinger, Nature Mat. 

11, 126 (2012). 

[26] A. Vanossi and E. Tosatti, Nature Mater. 11, 97 (2012). 
[27] A. Vanossi, N. Manini, E. Tosatti, Proc. Natl. Acad. Sci. 
U.S.A. 109, 16429 (2012). 



